Low-frequency characterization of quantum tunneling in flux qubits 
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We propose to investigate flux qubits by the impedance measurement technique (IMT), currently 
used to determine the current-phase relation in Josephson junctions. We analyze in detail the case 
of a high-quality tank circuit coupled to a persistent-current qubit, to which IMT was successfully 
applied in the classical regime. It is shown that low-frequency IMT can give considerable information 
about the level anticrossing, in particular the value of the tunneling amplitude. An interesting 
difference exists between applying the ac bias directly to the tank and indirectly via the qubit. In 
the latter case, a convenient way to find the degeneracy point in situ is described. Our design only 
involves existing technology, and its noise tolerance is quantitatively estimated to be realistic. 

PACS numbers: 74.50.+r, 84.37.+q, 03.67.-a 



I. INTRODUCTION 



Josephson-junction flux qubits are known to be^andi- 
dates for solid-state quantum computing circuits!!! This 
qubit variety has good tolerance to external noise, es=. 
pecially to dangerous background-charge fluctuations.tl 
A flux qubit is a superconducting loop, the two lowest- 
energy states of which differ in the direction of circulating 
persistent current. For many flux qubits, these two states 
become degenerate when the external flux <I> X threading 
the loop equals <E>o/2 ($o = h/2e is the flux quantum), 
and quantum tunneling between them becomes possible. 
Moving $ x away from $0/2 lifts the degeneracy and ap- 
plies a bias between the two states. When the biasing 
energy exceeds the tunneling amplitude A the tunneling 
stops, but the relative phase between the two states will 
still evolve in time. This, together with coherent tun- 
neling, provides single-bit quantum gate operations. To 
have a universal set of gates, necessary for quantum com- 
puting, one needs to be able to couple two qubits. The 
methods of coupling two flux qubits and performing gate 
operations are beyond the present scope. Instead, we pro- 
pose a method to characterize the quantum behavior of 
a flux qubit by coupling it to a tank circuit. The discus- 
sion will be quite general and can be applied to different 
types of flux qufrit such as rf-SQUIDjj three- Josephson- 
junction (3JJ),u'u multi-terminal,tj etc. We will use the 
example of the 3JJ qubit, where quantum superposition 
of the macroscopic current states has been observed^ 

Due to the loop self-inductance, the total qubit flux $ 
may differ from $ x , depending on the direction of the 
persistent current. Figure ^ shows the $-$ x curve for 
a typical flux qubit. The solid lines correspond to clas- 
sical behavior. Near the degeneracy point, the diagram 
is hysteretic, a signature of the qubit's bistability. This 
has been observed for the 3JJ in Refs. [?|,||. In the quan- 
tum regime, tunneling between the states at degeneracy 
may eliminate the hysteresis (dashed line in Fig. [j]) . This 
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FIG. 1: Solid lines: hysteretic dependence of the total flux $ 
on the external flux <I> X in the classical regime. Dashed line: 
disappearance of the hysteresis by quantum tunneling. 



phenomenon will be discussed in detail below. 

In general, one can plot the classical (local) minimum 
energies of a flux qubit as in Fig. |^a. The left (right) 
branch then corresponds to (counter-)clockwise flow of 
the spontaneous current. The hysteresis is also evident 
from this diagram. In the quantum regime, there will be 
discrete local states in each of the qubit's bistable po- 
tential wells. From now on we denote the lowest-lying 
such states as "J 1 and 5 ,r , corresponding to "left" and 
"right" directions of the persistent current respectively. 
At <& x = $0/2, resonant tunneling will render the low- 
est eigenstates of the full Hamiltonian as superpositions 
(^'i* 1 ') / 'V2. A small splitting equal to 2A will appear 
between their energies (Fig. ||b) . Starting with the qubit 
in its ground state (lower band in Fig. ^b), adiabatically 
changing <i> x will keep it in the ground state. This means 
that by passing through the degeneracy point, the qubit 
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FIG. 2: (a) Minimum energies of a qubit as a function of 
external magnetic flux in the classical regime, (b) Quantum 
mechanical energy profile for the same qubit as in (a). 



will continuously transform from 'F 1 to , F r . This pure 
quantum behavior is shown by the dashed line in Fig. [I]. 
On the other hand, if $ x changes rapidly, there is a con- 
siderable probability to excite the qubit and therefore 
continue on the same classical branch (left or right). This 
so-called Landau-Zener effect can be used to distinguish 
the classical from the quantum energy curves. 

The curvature of the energy profile is related to the 
qubit 's effective inductance and is therefore important 
for measurement. Figure || displays the second derivative 
of the curves in Fig. In the classical regime (Fig. [||a), 
the hysteretic behavior is the same as for the energy. On 
the other hand, in the quantum regime the hysteresis is 
replaced by a sharp spike, due to the level anticrossing. 
The appearance of this spike can be ascribed to enhanced 
susceptibility of the system due to tunneling. Its size and 
width can provide information about A. 



(a) 



FIG. 4: Flux qubit coupled to a tank: direct biasing scheme. 
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FIG. 3: (a) Second derivative of the qubit's classical minimum 
energy vs external magnetic flux, (b) Second derivative of the 
same qubit's ground-state energy in the quantum regime. 

A simple experimental implementation is to induc- 
tively couple the qubit to an LC tank circuit with known 
inductance Lt, capacitance Ct, and quality Q through 
a mutual inductance M (Fig. ^). The resonant charac- 
teristics of the tank circuit (frequency, phase shift, etc.) 
will then be sensitive to the qubit inductance and there- 
fore to its energy curvature. In particular, the spike in 
Fig. Bp appears as sharp dips in both phase shift and 



tank voltage as a function of <& x (see Section III ) 



This method, known as impedance measurement tech- 
nique (IMT), has been used for current -phase measure- 
ments of Josephson junctions. It originates from the pi- 
oneering work of Rifkin and Deaverp and is analyzed in 
detail in Ref. llOl IMT has also successfully been applied 



to a 3JJ qubit in the classical regime.0 and the hysteretic 
dependence of the ground-state energy on <I> X (cf. Fig. ||a) 
was observed as predicted in Ref. ^. The method has also 
been used for the investigation of quantum transitions in 
an rf-SQUID (Ref. [ll] and references therein). 

First of all, in Section || we calculate the two qubit 



energies in more detail than in Ref. m In Section III 



study the qubit's interaction with a high-quality resonant 
tank, showing that low-frequency IMT yields useful in- 
formation abo ut the qubit's quantum behavior. Finally, 
in Section IV , the effect of noise is considered. 



II. QUANTUM DYNAMICS OF THE 3JJ QUBIT 

The 3JJ qubitfi consists of three Josephson junctions in 
a loop with very small inductance L, typically in the pH 
range. This insures effective decoupling from the environ- 
ment. Two junctions have equal critical current I c and 
(effective) capacitance C, while those of the third junc- 
tion are slightly smaller: al c and aC, with 0.5 < a < 1. 
If the Josephson energy Ej = I c $>q/2tt is much larger 
than the Coulomb energy Ec — e 2 /2C, the Josephson 
phase is well defined. Near $ x| -jp-.$o/2, this system has 
two low-lying quantum states.u!13 The energy splitting 
between them in the presence of a small flux bias has 
been given in Ref. q but only for a particular choice of, 
e.g., a and g = Ej/Ec- In this section we derive the 
splitting with its explicit dependence on the qubit pa- 
rameters. The energy levels are derived from the Hamil- 
tonian [see Eq. (12) in Ref. [| 
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where ip = (ipi+(p 2 )/2, 9 = (ipi~ip 2 )/2 with ipi t2 
the phase differences across the two identical junctions, 
ihd v , Pg = -ihde, M v = ($o/2vr) 2 2C, M e = 
and 



P v = 



(l+2a)M v , 

U(f x ,ip,6) = Ej[a-2costpcos6 + acos(2irf x +29)} . (2) 

In contrast to Ref. we define the flux bias / x = 
$x/$o — 5 as a small deviation from degeneracy. 

Since the qubit is assumed to have small L and J c 
(typically L w 10 pH, I c w 100 nA), the shielding factor 
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LI c /$q w 0.001. Hence, in (|]) we have neglected the 
shielding current, considering <I> as an external flux. 

At / x = 0, the potential (^) has two minima at (p = 0, 
= ±0», with cos6>* = l/2a (0, > 0). Tunneling lifts 
their degeneracy, leading to energy levels E± — £o ± A. 
To find the levels for | / x | <C 1 we expand Eq. (||) near its 
minima, retaining linear terms in / x and quadratic terms 

r/l 

in 4>, 9. Define 6 J as the minima, shifted due to / x : 



,r/i 



= ±9, + 27r/ x 
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(3) 



that is, the upper (lower) sign refers to the right (left) 
well. The potential energy then reads 



U 
Eh 



2a a 2a \ V4a 2 -1 



Q r /h 2 



2a-±±f^^L 
2a a V4a 2 -1 



(4) 



Near degeneracy, the eigensolutions of Hq^± = E± 1 i± 
can be written as superpositions ^± = a±^ 1 + b±^> T , 
yielding the well-known eigenenergies E± = (e'+e r )/2 ± 
v /(e 1 -e 1 ') 2 /4 + A 2 , with e r /' = (W^\H \^ t/1 ). The ma- 
trix element A cannot accurately be found in terms 
of W/ x . In what follows it is assumed constant, 
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(5) 



neglecting its dependence on / x . 

To find the dependence of E± on / x , we take \E' r / 1 to 
be oscillator ground states in their respective wells: 
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V «5 V (4a 2 -l) 3/2 , 



Combining the above, one finds the eigenenergies 



E ± =e ±jE 2 fi\i(a) + A2 



(6) 

(7) 

(8) 
(9) 

(10) 



where 
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V4a 2 -1 



(11) 



(12) 



The splitting given by Eq. ( |lO| ) differs from that of Eq. (1) 
in Ref. || by a factor A (a) which explicitly accounts for 
the dependence of E± on a and g. 

For stationary states, the current in the qubit loop 
can be calculated either as the average of the current 
operator I q — I c sm(ip + 9) over the eigenf unctions, or as 
the derivative of the energy over the external flux: 

I q = (* ± |/ q |^) = ^=±/ c / x ^^, (13) 
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where hu>o = E+ — E- . In equilibrium at finite temper- 
ature T, Eq. (H3) readily generalizes to 



i q = (* + |/ q |*+>/>?V + (*_|/ q |*_> P !?_ 



2k B T 



(14) 



with the density matrix elements p+ + = e E +/ ksT /Z and 

p aq_ = e -E-/k B T /z ^ where Z = e -E+/k B T + e -E-/ kB T_ 



III. QUBIT-TANK INTERACTION 

We propose here to extract information about the 
quantum dynamics of a flux qubit with the aid of a clas- 
sical linear high-quality tank circuit, coupled to the qubit 
via a mutual inductance M. The tank consists of a ca- 
pacitor Ct, inductor Lt, and a resistor Rt which are 
connected in parallel and driven by a current source Jb (t) 
(Fig. ||). The problem of coupling a quantum object to 
a dissipative classical one has no unique theoretical so- 
lution. However, if we assume that the classical object 
is much slower than the quantum one we may solve for 
the latter's motion, accounting for the coupling-coordi- 
nates of the former as mere external parameters.liiJ Here, 
the characteristic frequency A/ft of the qubit is in the 
GHz range, while the resonances wt of our tank circuit 
lie below 100 MHz. There exist two different schemes of 
coupling a tank circuit to the qubit. First we consider 
direct biasing, where a current Ib{t) = Iq cos tot is fed 
directly into L T (Fig. 0). 



A. Direct biasing scheme 



The voltage across the tank circuit evolves as 

V + ^-V + uo T V = -Mu T i q + -^i b {t) . 



(15) 
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FIG. 5: Tank phase x vs bias amplitude Jo; A/ft, = 300 MHz. FIG. 6: Tank phase \ vs bias flux / x ; Iq = 100 pA. From the 
From the lower to the upper curve, the bias flux 10 4 / x takes lower to the upper curve (at / x = 0) , the tunneling frequency 
the values 0, 2, 4, 6, 8, 10. A/h takes the values 150, 300, 450, 600, 750 MHz. 



Here, Q = lotRtCt ^> 1 and lot — I/V-^tCt! Iq is 
given by ( |i~3| ) or (|l4|), and depends on the qubit flux 
$ = $ x + MI L , where I L = L^ 1 J V dt is the current 
in Lt and $ x is time independent. Below we study the 
simplest case k^/T <C A, so that the qubit is definitely in 
its ground state E_. Then, Eq. (J15|) takes the form 



V 



^-V + u 2 V 



where k 2 = M 2 /LL T , 
d 2 E- 



-k 2 Lu 2 T ^-V 
d$ 2 



I 2 A 2 X 2 (a 
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(2tt) 2 (£Jf/ 2 A 2 (a) + A 2 ) 



3/2 



j 6 (t) (16) 



(17) 



and / = [$ x +M/i(t)]/$ - | (1/1 < 1)- Thus, © 
is nonlinear in V. Since the coupling to the qubit is 
small one may apply the method of hamionic balance, 
which is well known in rf-SQUID theory.E3 Accordingly, 
if uj rs wt, then V oscillates with frequency lj, while 
its amplitude v and phase x are slow functions of time: 
V(t) = v(t) cos[ujt + x{t)]- From (111) we obtain 



X = ^tCo 
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with the detuning £o = (wt _ w)/wt, and where 



1 f 2 " 

i p («./x) = - / d4>- 

n Jo \i + n 



cos" 



[l + V 2 (/x + 7 sin </>)'] 



2l 3/2 



, (20) 



with 77 = _EjA(ev)/A and 7 = Mv/ojtLt'&o- 

Setting u = x = in (|18|) and ( |i~9| ) one obtains the 
stationary tank voltage and phase, 

w 2 (1 + 4Q 2 £>, / x )) = /o 2 ^ 2 L^Q 2 (21) 
tan X = 2Q£( V ,/x) , (22) 
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FIG. 7: Tank voltage v vs bias flux / x ; I = 100 pA. From the 
lower to the upper curve (at / x = 0), the tunneling frequency 
A/h takes the values 150, 300, 450, 600, 750 MHz. 



where we introduced a flux-dependent detuning 



LI 2 fX(a) 



2tt 



(23) 



We have used Eqs. (|2l|)-(|23|) to find voltage-flux v(f x ), 
phase-current x(^o), and phase-flux x(/x) characteristics 
at resonance uj = oj^. We take the qubit parameters 
I c = 400 nA, a = 0.8, L = 15 pH, g = 100, a tank 
with L T = 50 nH, Q = 1000, lu t /2tt = 30 MHz, and 
k = 10 -2 . The x(-^o) curves for several / x are shown in 
Fig. |. The x(/x) and v(f x ) curves are shown in Figs. 
U and ^ for various A. The sharp dips in Figs, [fj] and 
correspond to the spike in the second derivative of the 
energy profile in Fig. ||b. Clearly, the width of the dips 
is correlated with A: with the increase of A the width 
of the dips also increases. The x(/x) curves for different 
Io are shown in Fig. ||. The shape and the value of x are 
seen to be very sensitive to Iq. The dependence of the 
voltage modulation 6v = v(f x —0) — v(f x —10^ 3 ) on Iq is 
shown in Fig. |^. 



5 




-0.003 -0.002 -0.001 0.000 0.001 0.002 0.003 

FIG. 8: Tank phase X vs bias flux / x ; A/h = 150 MHz. From 
the lower to the upper curve (at / x = 0), the bias amplitude 
Jo takes the values 10, 50, 100, 150, 250 pA. 
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FIG. 9: Voltage modulation 5v = w(/ x =0) - u(/ x =10~ 3 ) vs 
bias current 7 ; A/A = 300 MHz. 



B. Scheme with separate driving coil 

In this scheme, a bias flux $b (i) = 'J'ac sin o>i is applied 
to the qubit loop from a separate coil (Fig. |h]). The tank 
response is similar to (pi): 



V+^V + uj^V = 



1 



-Mli) T — $ac + ^ac WCOSO>£, 

(24) 

where $ a c is the flux which the external coil couples 
directly into the tank and d^,E_ is given by (17) with 
/ = ($ x +$ac sin^)/$ - | = fx + /ac sin o>t. Rewriting 



the first term on the rhs of (24) as 



-k^y/I^L c4w* /acG(t) (25) 



makes its time dependence manifest: 
G(t) = C08wt 



[l + rf (f x + f ac smu)t) 



2-1 3/2 



(26) 



The advantage of a separate driving coil is that one 
can effectively decouple the tank from the fundamental 
harmonic of the bias, since the qubit signal G(t) contains 
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FIG. 10: Flux qubit coupled to a tank: scheme with a separate 
driving coil. 
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FIG. 11: Tank voltage vs bias flux / x ; A/h = 300 MHz, 
"I'ac = 5 • 10~ 4( I>o. Dotted line: driving frequency lu = lut/2; 
solid line: u> = wt/3. 



not only uj but also 2u, 3a;, etc. At / x = 0, G(t) contains 
only odd harmonics. This can be used to find the degen- 
eracy point in practice, e.g., by tuning / x so that the tank 
response vanishes (reaches its maximum) at frequency 2o> 
(3a;). We have studied the higher harmonics by solving 
Eq. (|24| ) numerically with <f> ac = for to — ujt/2 and 
uj = wt/3, see Fig. Since the full amplitudes contain 
contributions from all harmonics, at u> = o>t/2, / x = 
one observes a finite dip rather than a zero. 



IV. REQUIREMENTS ON NOISE SOURCES 

Figures || and clearly reveal the quantum nature of 
the flux qubit within a range \Sf x \ < 5 • 1CP 4 from the 
degeneracy point / x = 0. Therefore, the unavoidable 
external flux noise coupled to the qubit must be much 
smaller than this value. The most i mportant sources are 
the Nyquist noise I n — ^ 'Ak^T/ Rt and the current noise 
7 a of the preamplifier. The former generates the qubit- 
flux noise <I> n = MI n Q\fB, where B = ujt/2-kQ is the 
tank bandwidth. With T = 20 mK and the tank pa- 

- 6 $n 



rameters of Section III, one gets $ n s» 8 • 10 F° r 
7 a = 10 _14 A/\/Hz, we estimate the corresponding flux 
noise as <f> a = MI a Q\/~B w 7 • 10 _6 $o- Thus, the noise 
these sources couple to the qubit is at least two orders 
smaller than the peak widths in Figs. ^ and [?]. On the 
other hand, these sources give rise to directly detected 
voltage noise across the tank circuit. The thermal tank 
noise is V n = I^uj^L^Q^/B w 17.6 nV. The noise due 
to 7 a is Va,i — IhUJtLtQVB w 16 nV. And finally, if we 
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take Va.2 — 40 pV/vHz for the preamplifier's own volt- 
age noise, we get V^VB ss 7 nV for its contribution in 
the tank bandwidth. Comparing these values with the 
voltage modulation in Figs. [?], || and [ll], we see they are 
at least several times smaller than the qubit signal. 

V. CONCLUSION 

We have shown that IMT can be used for low-frequency 
characterization of the ground (in general: equilibrium) 
state of a flux qubit. The method allows determining the 
tunnel splitting between qubit states for a broad class 
of devices; with the term "flux (as opposed to phase) 
qubit" we stress that the two states must differ not only 



in Josephson phase, but in actual magnetic flux visible 
to the outside] 1 ] The design exclusively employs present- 
day technology, and the expected noise levels have been 
shown not to disrupt the measurement. On the qubit 
time scale, the method is a quasi-equilibrium one; as 
such,Ji can determine A but not, e.g., the "T2" dephasing 
time.li-3 Efforts to adapt IMT to this and related problems 
such as qubit readout and control are underway. 
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